Errors in Estimating Q\ due to the Fluid Approximation 
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The matter content of the Universe is strongly inhomogeneous on small scales. Motivated by 
this fact, we consider a model of the Universe that has regularly spaced discrete masses, rather 
than a continuous fluid. The optical properties of such space-times can differ considerably from the 
continuous fluid case, even if the 'average' dynamics are the same. We show that these differences 
have consequences for cosmological parameter estimation, and that fitting to recent supernovae 
observations gives a correction to the inferred value of S7a of ~ 10%. 



In the standard approach to relativistic cosmology it is 
usual to assume that the cosmological principle implies 
that the Universe is permeated by a set of continuous 
perfect fluids. Examples of these are photons, baryons, 
dark matter and dark energy. Yet when we observe the 
Universe, there is clear evidence for discreteness. Mat- 
ter is accumulated in stars and galaxies which arrange 
themselves in clusters, filaments and walls. These struc- 
tures occupy a small volume of space, the rest of which 
is almost completely devoid of electromagnetically inter- 
acting matter. It is therefore reasonable to ask if the 
fluid approximation is a good representation of the real 
Universe, and, if not, what are the corrections due to the 
discretization we observe. 

In [l[ (henceforth, paper I) we considered just such a 
model of the Universe: The matter content was taken to 
be discrete islands of mass, rather than the usual contin- 
uous fluid. This approach, which we dubbed 'Archipela- 
gian Cosmology', was based on the Lattice Universe 
model of Lindquist and Wheeler 0] , a construction anal- 
ogous to that of Wigner and Seitz in electromagnetism 
[3| . The principal result of the study was that even if the 
large-scale dynamics of a universe with discretized mat- 
ter content approach those of a Friedmann-Robertson- 
Walker (FRW) universe, this does not mean that the op- 
tical properties will. Observers in a universe filled with 
discrete objects can measure different redshifts and lumi- 
nosity distances to distant astrophysical objects, even if 
the 'average' expansion is identical to FRW. 

Of course, paper I is not the first study of the optical 
properties of an inhomogeneous universe. Previous at- 
tempts have been made in the context of Swiss Cheese 
cosmologies 0, H, 0, H, H , as well as with linear per- 
turbations to FRW geometry [H, El El, El, El- The 
difference in paper I is that the model used does not rely 
on an FRW background (neither as an embedding struc- 
ture for the inhomogeneities, nor as a background to be 
perturbed around). Instead, the model is a bottom up 
construction that is not FRW at any point in space-time, 
yet behaves dynamically like FRW on large scales. The 
model does not, therefore, rely on any ill-defined concepts 
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such as 'average' geometry, and is free from ambiguities 
associated with the scale of smoothing. Hence, we ex- 
pect this approach to provide us with new insights, and, 
in particular, new ways to test concepts such as the fluid 
approximation without the constraint of being tied to 
FRW geometry. 

In paper I we considered the simple case of Milky Way- 
sized spherically symmetric masses arranged on a regu- 
lar lattice with a spacing of ~ IMpc, and with critical 
density. We found that in the absence of a Cosmolog- 
ical Constant the deceleration parameter measured by 
an observer in such a space-time would be qa ~ 8/7, 
rather than the usual value of 1/2 that would be mea- 
sured in a universe filled with a perfect fluid of dust (i.e. 
an Einstcin-de Sitter universe). This result is due to two 
separate effects. Firstly, the null geodesies that connect 
observers and sources pass through the empty regions be- 
tween masses, rather than through a continuous energy 
density. This means that the focusing due to the notional 
fluid of an FRW cosmology is absent [H, EH. Secondly, 
the redshift experienced by the photons is due to the 
effect of anisotropic expansion integrated along their tra- 
jectories. This does not, in general, correspond to the 
same redshift that would be experienced by a photon in 
the 'average' space-time. Analytic expressions for these 
quantities are given in paper I, together with the results 
of numerical simulations. 

In this paper we wish to include the effect of a Cosmo- 
logical Constant. In some sense, this is a very straightfor- 
ward extension of the cosmology that was considered in 
paper I. Crucially, however, it will allow us to assess the 
influence of the effects we have uncovered on cosmological 
parameter estimation, as well as provide us with a test 
of the validity of our approach. Huge resources are be- 
ing invested into gathering ever more observational data, 
and it is critical to the success of the missions involved in 
this that we understand to the highest degree possible the 
relationships between redshifts, luminosity distances and 
the expansion of the Universe. Any attempts at precision 
parameter estimation [l7l |. or, for example, reconstruct- 
ing the equation of state of Dark Energy [l8j], will be 
highly dependent on such considerations. 

Let us consider a critical density cosmology where we 
have a lattice of cells with Schwarzschild-de Sitter ge- 
ometry. The line-element inside each cell is then given 
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Here the time coordinate r is the proper time of a freely 
falling observer with 4-velocity 
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The relationship between r and the Schwarzschild time 
coordinate t is given in Appendix A of paper I. The reader 
is also referred to Section 3 of paper I for an explanation 
of why t corresponds to 'cosmological time'. 

Now, the principal difference between the Lindquist- 
Wheeler lattice construction in General Relativity (GR) 
2], and Wigner-Seitz construction in electromagnetism 
3j (aside from the non-linearity of the field equations 
in GR) is that in the case of the former the lattice itself 
becomes dynamical. In the case of a universe with critical 
density, this leads to cosmological evolution governed by 
the equation 
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where the over-dot here denotes differentiation with re- 
spect to r. If we replace r by the scale factor a, then this 
is clearly just the Friedmann equation, with the usual 
solution 
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where L$ is the lattice spacing today (so that a(ro) = 1). 
The large-scale dynamics of this model are therefore iden- 
tical to an FRW model with 8np — 2m/ (Loa) 3 . However, 
as discussed above, this does not necessarily mean that 
the optical properties are also the same. 

We know from our study of the ^a = case that these 
corrections can be considerable. However, we also know 
that as f^A — > 1 they should vanish. This is due to the 
space-time approaching de Sitter space in this limit. The 
influence of the mass at the centre of each cell then be- 
comes negligible, as the space-time is dominated by A. Of 
course, the energy density associated with A is constant 
everywhere, and so we should expect any effects due to 
the discreteness of m to disappear. Indeed, we find just 
such a convergence to occur in the results below, and we 
consider this to be an important verification of the model 
we are using (analogous to Shockley's 'empty lattice test' 
with a constant potential fl9j |, that is used to justify the 
Wigner-Seitz approximation) . 

Now, to begin our study of redshifts and distance mea- 
sures in this space-time, let us consider a bundle of null 



rays that connect an observer with an unobscured source. 
The null geodesic equations can then be integrated to give 
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where B, J and are constants, and over-dots here 
correspond to derivatives with respect to the affine pa- 
rameter, A. These equations allow us to propagate null 
geodesies away from our observer, and can be connected 
between lattice cells using the methods discussed in pa- 
per I (with 2m j r replaced by 2m jr + Ar 2 /3, wherever it 
occurs in the matching conditions). 

The redshift between emission and observation is 
given, as always, by the ratio of frequencies at these two 
events. For source and observer both moving with 4- 
velocity u a , as given in ©, this is 
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where k a is the 4- vector tangent to the null geodesic, 
and subscripts e and o denote quantities at emission and 
observation, respectively. The null geodesies equations, 
and redshift, can now be calculated numerically. 

For the = case we obtained in paper I the analytic 
approximation 



l + z~ (l + z F w) <7> , 
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where (7) = 0.7. However, we know that in the limit 
ft A — * 1 that we should find (7) — > 1, as the space-time 
approaches de Sitter space. We therefore expect that we 
should find 0.7 < (7) < 1, for < fl A < 1. 

Our numerical simulations verify this expectation, and 
we find that excellent fits can be achieved in the range 
< z < 2 by allowing (7) to be a running function of 
zfrWi with the form 



(7) = A + Bzfrw- 



(11) 



Performing a least squares fit to the output of our numer- 
ical code we find best fit values of A and B for various 
different values of These are displayed in Fig. [IJ 
along with the fitted functions A = 0.69 + 0.29J7a and 
B = 0.0021-0.057O\ + 0.055f2 A 2 . It can be seen that (7) 
approaches a constant as both Q\ — > (as found in paper 
I), and as fl\ — > 1 (as the de Sitter limit is approached). 

To go further, and obtain expressions for angular di- 
ameter, and luminosity distance in these models, we 
must integrate the Sachs optical equations [2fJ along the 
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FIG. 1: Dots show the best fit values of A and B from Eq. 

for various values of Qa- The curves in the upper and 
lower plots are the fitted functions 0.69 + 0.29J1a and 0.0021 — 
0.057f2 A + 0.055f2 A 2 , respectively. 

geodesies that are solutions of ([H])-©- Without rotation, 
these equations are 
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^ + 2a§ = C abcd (t*) a k b (t*yk d = C, (13) 
dX 



where 9 and a are the expansion and complex shear 
scalars, respectively. The C a bcd is Weyl's tensor, R a b 
is the Ricci tensor, and t a is a vector that is orthogonal 
to k a , null, and has a magnitude of 1 (i.e. t a k a = 0, 
t a t a — and t a (t*) a — 1). The initial conditions for 
integrating these equations are a\ Q = 0, ta\o = and 
drA/d\\ — constant. 

Once the expansion scalar is known, then the angular 
diameter distance is given by the integral 



ta oc exp 
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and the luminosity distance is given by Etherington's the- 



orem [21| as Th = (1 



ta- For the geometry specified 



by (jTJ) we then find that the driving terms in (fT5|) and 
(fTBf are given by 
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where "J is a constant, specifying the complex phase. 
These equations are independent of Q\, and so their so- 
lution is very similar to the fl\ — case (up to the effect 
JIa has on the trajectories r = r(X)). 

By numerically integrating (fT2|) and (fl~3|) we find that 
the shear is typically unimportant at z £ 1, unless a 
trajectory happens to pass very close to a central mass. 
At higher redshifts, however, the shear can accumulate 
and become more important, eventually causing the di- 
vergence of the optical scalars along some trajectories. 
This corresponds to caustics occurring in the beams. 



For the present purposes, at z < 1, it is sufficient to 
neglect the shear, and set a = 0. The solution to Eqs. 
(fT2l) and (fT4|) is then given simply by ta oc A, and so 
we find that rj, oc (1 + z) 2 X. We now need to find the 
relationship between A and z. Using the expression for 
redshift, (|10|) . together with the Friedmann equation, J3J, 
we can write 
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where the usual expressions for fi m and Oa have been 
used. In general, this equation needs to be integrated 
numerically to find a solution for X(z). However, we find 
that by treating (7) as a quasi-static variable we can 
obtain analytic results that are accurate to better than 
1%. In this case, the solution to Eq. (|T6| is 
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and 2-Fi(a, b; c; d) is the hypergeometric function [25j. We 
then have that, in the absence of shear, the luminosity 
distance is well approximated as a function of z by 



r L cx (l + z) 2 



L(i 



3 + 2( 7 ) 
z) 2 <"<> 



/(0) 



(17) 



We will now fit our model to the recent supernova ob- 
servations, and obtain an estimate for the difference in 
the best fit value of £Ia due to using a discretized matter 
content, rather than a perfect fluid. This will be done 
using the first-year Supernovae Le gacy Survey (SNLS) 
data, consisting of 115 supernovae [22j, calibrated with 
the Spectral Adaptive Light-curve Template (SALT) fit- 
ter [23| . Assuming a critical density, both FRW and 
our model have 5 free parameters: £Ia and the 4 'nui- 
sance' parameters required to calibrate the supernova 
data. These are the absolute magnitude, the intrinsic 
error and the color and stretch parameters used in the 
process of light curve fitting, {Mo, a int ,a, (3}. 

We show the best fits for both our model, and a spa- 
tially flat continuous fluid FRW cosmology in Fig. [2 The 
fits are very similar with a difference in log likelihood of 
just Aln£ = 0.37, in favor of the FRW model. In fact, 
it can be seen from the figure that the distance moduli of 
these two best fit models are virtually identical, with the 
two solid lines only actually becoming distinguishable by 
eye at z > 0.8. It seems clear that a lot more data will be 
needed before there is any hope of distinguishing these 
two curves observationally. 
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FIG. 2: The best fit critical density models with discrete matter (thick line) and continuous fluid (thin line) sources. The data 
is the 115 high and low redshift supernovae from SNLS [22|. fitted to the discrete matter model with the SALT light curve fitter 
[23l ]. Einstein-de Sitter (dotted line) and de Sitter (dashed line) models are also shown, for reference. The distance modulus, 
Adm, is the magnitude of the source, minus the magnitude it would have at the same 2 in an empty, open Milne universe. 



The interesting point, and the main result of this pa- 
per, is that these two best fit models correspond to 
quite different values of Qa- For the more usual per- 
fect fluid FRW cosmology we recover the standard result 
qFRW _ 0.74 ± 0.04. For the model with a matter con- 
tent composed of discrete masses, however, we obtain 
SIa = 0.66 ± 0.04. The best fit values of the nuisance 
parameters are very similar in each would be 

expected from the similarity of the two distance moduli. 

One may initially react by observing that this differ- 
ence in £7a is only at the level of 2a, and as such is not 
very significant. Our point, however, is that irrespective 
of how good the data eventually becomes, there will still 
be a difference in Qa of ~ 10% between these two mod- 
els. This difference is not due to any new physics, exotic 
matter content, or unexpectedly large structures of voids 
in the universe, but only to the fact that we have treated 
the matter content of the Universe in our model as being 
discrete, rather than continuous. 

Inevitably the reader will be concerned with the gen- 
erality of this result, and whether it should apply to the 
real Universe, or only to our simple model. We strongly 
suspect that the real Universe, with its complicated net- 
works of voids, filaments, walls and nodes [24[ will not 
behave exactly as our model does, and that there will 
be new behaviour beyond that which we have uncovered 
here. In future work we hope to make progress in gener- 
alizing our results to more realistic situations. However, 
we believe that there is good reason to expect some of 
the essential features of our study to hold in the more 
general case. 

Firstly, the photons we observe on Earth have not 
reached us by travelling through a continuous fluid 
of critical density; they have travelled through mostly 
empty space. As such, they should not experience the 
focusing that such a fluid would produce. Photons ex- 
perience the geometry of space-time through which they 
have passed, and not the global average. Secondly, the 
expansion that photons experience is also not the global 



average, it is the integrated effect of the expansion which 
is strongly anisotropic at any given point along their tra- 
jectory (as long as il m > 0). For further discussion of 
these points, and discussion on the appropriateness of 
using an approximate solution of Einstein's equations, 
we refer the reader to paper I. 

Of course, one will also be concerned about system- 
atic errors, and the validity of the assumptions that have 
gone into this model. In future publications we intend to 
establish the validity of these assumptions through de- 
tailed investigation, but for now we already good reason 
to suspect that the approximations that have been made 
are good. In paper I it was shown that the approximate 
nature of the boundary conditions between lattice cells 
is insensitive to the details of the matching conditions, 
lending credence to the idea of an average tangency be- 
tween neighbouring space-like surfaces. Furthermore, we 
also obtained in paper I analytic approximations that are 
in good agreement with our numerical results. These al- 
low some physical insight into the source of the effects 
we have uncovered, and further their legitimacy. Beyond 
this, the analogy with the highly successful Wigner-Seitz 
construction and the recovery of the familiar optics of de 
Sitter space in the appropriate limit is also encouraging. 

In summary, we find that the fluid approximation 
in cosmology, while appearing innocuous, can introduce 
considerable errors in interpreting cosmological data. Us- 
ing a simple model of the Universe, with discrete masses 
arranged on a regular lattice, we have shown that even if 
the average dynamics of the Universe are unchanged, the 
error introduced in the estimation of fl\ due to different 
optical properties can be of the order of 10%. Such ef- 
fects will need to be understood and accounted for if we 
are to attempt precision cosmology. 

Acknowledgements 

We are very grateful to George Ellis, Chris Clarkson 
and Lance Miller for suggestions and discussion. TC ac- 
knowledges the support of Jesus College and the BIPAC. 



5 



[1] Clifton, T. & Ferreira, P. G., Phys. Rev. D80, 103503 
(2009). 

[2] Lindquist, R. W. & Wheeler, J. A., Rev. Mod. Phys. 29, 
432 (1957); erratum, Rev. Mod. Phys. 31, 839 (1959). 

[3] Wigner, E. P. & Seitz, F., Phys . Rev. 43, 804 (1933). 
Wigner, E. P. & Seitz, F., Phys. Rev. 46, 509 (1934). 

[4] Kantowski, R., Astrophys. J. 155, 89 (1969). 

[5] Biswas, T., Mansouri, R. & Notari, A., J CAP 12, 17 
(2007). 

[6] Biswas, T. & Notari, A., JCAP 6, 21 (2008). 

[7] Brouzakis, N., Tetradis, N. & Tzavara E., JCAP 2, 13 

(2007); JCAP A, 8 (2008). 
[8] Marra, V., Kolb, E. W., Matarrese, S. & Riotto, A., 

Phys. Rev. D76, 123004 (2007). 

[9] Clifton, T. & Zuntz, J., |arXiv:0902.0726] (2009). 
[10] Futamase, T. & Sasaki, M., Phys. Rev. D 40, 2502 

(1989). 

[11] Watanabe, K. & Tomita, K., Astrophys. J. 355, 1 (1990). 
[12] Kasai, M., Futamese, T. & Takahara, F., Phys. Lett. A 
147, 97 (1990). 



[13] Holz, D. E. & Wald, R. M., Phys. Rev. D58, 063501 
(1998). 

[14] Kainulainen, K. & Marra, V., arXiv:0909.0822] (2009). 

[15] Bertotti, B., Proc. Roy. Soc. Lond. A 294, 195 (1966). 

[16] Dyer, C. C. & Roeder, R. C, Astrophys. J. 174, L115 
(1972); Astrophys. J. 180, L31 (1973). 

[17] Bridle, S. L., Lahav, O. & Ostriker, J. P., Science 299, 
1532 (2003). 

[18] Sahni, V. & Starobinsky, A., Int. J. Mod. Phys. D 15, 

2105 (2006). 
[19] Shockley, H., Phys. Rev. 52, 286 (1937). 
[20] Sachs, R., Proc. Roy. Soc. Lond. A 264, 309 (1961). 
[21] Etherington, I. M. H., Phil. Mag. ser. 715, 761 (1933). 
[22] Astier, P. et ai, Astron. Astrophys. 447, 31 (2006). 
[23] Guy, J. et al, Astron. Astrophys. 443, 781 (2005). 
[24] Sylos Labini, F., Vasilyev, N. & Baryshev, Y., A & A 

496, 7 (2009). 

[25] The solution of z{l-z)f"(z)+[c-{a+b+l)z]f'{z) =abf. 



